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Abstract 

In any spatially discrete model of pedestrian motion which uses a regular 
lattice as basis, there is the question of how the symmetry between the dif- 
ferent directions of motion can be restored as far as possible but with limited 
computational effort. This question is equivalent to the question ''How im- 
portant is the orientation of the axis of discretization for the result of the sim- 
ulation ? " An optimization in terms of symmetry can be combined with the 
implementation of higher and heterogeniously distributed walking speeds by 
representing different walking speeds via different amounts of cells an agent 
may move during one round. Therefore all different possible neighborhoods 
for speeds up to v = 10 (cells per round) will be examined for the amount 
of deviation from radial symmetry. Simple criteria will be stated which will 
allow to find an optimal neighborhood for each speed. It will be shown that 
following these criteria even the best mixture of steps in Moore and von Neu- 
mann neighborhoods is unable to reproduce the optimal neighborhood for a 
speed as low as 4. 
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1 Introduction 




Figure 1: von Neumann- and Moore-neighborhood (v=l) 



In a model which is spatially and temporally discrete the speed of an agent (as 
the model of a person in a simulation will be called) is the number of cells which 
he moves during one round. As the real- world interpretation of the size of a cell 
is fixed by the scale of the discretization, the real-time interpretation of one round 
fixes the real- world interpretation of such a dimensionless speed. 
Subsequent steps within one of the neighborhoods of Fig. [H leave an agent to be 
either \/2 times as slow or as fast moving into the diagonal direction than moving 
horizontally or vertically. To some extent the situation can be improved by a mix- 
ture of von Neumann and Moore neighborhoods. If for example an agent would be 
allowed to do a total of five steps during one round, of which three are in von Neu- 
mann and two in Moore neighborhood, a larger total neighborhood of cells which 
can be reached during one round would result. The question is: Is there an optimal 
total neighborhood for a given speed? And can it be composed of von Neumann 
and Moore neighborhoods? In vertical and horizontal direction there is no doubt 
about the neighborhood: For a speed v = Vm the neighborhood contains the cell 
of the agent and Vm cells in each horizontal or vertical direction. For any other 
direction there are cells for which it is not obvious if they should be part of the 
neighborhood. At the very beginning f or i; = 1 there is the question if one should 
use the von Neumann or Moore neigborhood. (See Fig. [T]). 

Definition: Complete neighborhoods are four-fold symmetrical neighborhoods 
where all cells which belong to the neighborhood are closer to the center cell than 
those which do not. 

Example: There are three complete neighborhoods for each v = 2 and v = 3. (See 
Figs.[21and[3l) 
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Figure 2: Complete neighborhoods for v=2 



Obviously one can limit the search for an optimal neighborhood to complete neigh- 
borhoods. 

The question is, which complete neighborhood represents the corresponding speed 
(which is number of cells in horizontal and vertical direction) best. However there 
might exist other alternatives, the criteria chosen here are such that discretization 
effects concerning the axis of discretization of the original plan are minimized. 
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Figure 3: Complete neighborhoods for v=3 



Therefore at first for each complete neighborhood the speed v{(f)) into each direc- 
tion has to be written down. 
Then the direction-averaged speed is calculated: 






(1) 



After that the squared deviation of speeds into each direction from this average is 
calculated. 



Av^J^J{v{^)-Vav)^d^ 



(2) 



The criteria for an optimal neighborhood are then 

• The direction averaged speed should be close to the corresponding integer. 

• The deviation from this average into different directions should be small. 



Complete neighborhoods up to i; = 10: The neighborhoods are named by the 
maximum squared distance of a cell from the center, implying that neighborhood Y 
includes all cells of neighborhood X <Y. For reasons of clarity in the following 
table the numbers are only shown in the second octant. 
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NB: Of neighborhoods in the previous table only the neighborhoods 1, 2, 4, 5, 8, 
10, 13, 17, 20, 29, 34, 40, 45, 58, 80 and 97 can be composed of the corresponding 
number of subsequent steps within von Neumann and Moore neighborhoods. (Take 
N von Neumann and M — v — N Moore steps and check how the largest possible 
neighborhood looks like.) 



1.1 v((j)) - the variation of speed with the direction of motion 

Since all complete neighborhoods have a fourfold axe- symmetry it is sufficient to 
calculate i;(0) f or < < 7r/4. 

v{(l)) is continuously composed from different functions resulting from different 
ranges of 0. The structure of those ranges depends on the shape of the edge of the 
neighborhood (See Fig. ^. 
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Figure 4: Example for calculating v{(j)) for one of the v=3 neighborhoods 



Definitions: 



5xi horizontal distance (in cells) of a border cell to the origin (the 
black cell in Fig. |4]); with i starting with at = 
(see the numbering in the black circles in Fig. |4]) 
5yi vertical distance of a border cell to the origin; with i starting 

with at = 
D a large distance 
AX Dcos{^) 
AY Dsm{(f)) 
N total (minimal) number of steps to reach a cell in distance D 

into direction arctan(^) <(/){< arctan( ^^'^^ ) 
n number of steps into direction (0/0) -^ {Sxi/6yi) 
N — n number of steps into direction (0/0) -^ (fe^+i/5?/^+i) 
To reach the point (AX/ Ay) an agent has to do n times a {Sxi/6yi) step and 
{N — n) times a (fe^+i/^^^+i) step. 
Such that 

(3) 
(4) 



n 6xi + {N — n) fe^+i = 


AX 


n 6yi + (TV - n) Syi^i = 


Ay 


Solving this for N leads to 




Ay-rAX 

6yi+i - rSxi+i 




where 

_ 6yi^i - 6yi 





(5) 

oxi^i - dxi 

In the intervall [0, 7r/4] r can only take the values oc (range A in Fig. |4]) and -1 

(range B) which in the latter case in equation ([5]) has to be understood as limit, r is 

the local gradient of the border of the neighborhood. 

Since speed is distance (in cells) over number of rounds to move that distance, one 

has 

^ D ^ h.» - rfa.^. 

N sm (f) — r cos cp 
See Fig. [5]for the speed's dependence on the direction of motion of all three com- 
plete V = 2 neighborhoods. 



1.2 The integrals for the average speeds and the deviations 

For ranges with the same gradient of the border as in range A (vertical, see Fig. ^ 
to get the average one has to integrate: 

If = / -# = In ^ ^ (8) 



Speed dependence on direction of motion 
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for v=2 neighbourhoods 
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Figure 5: Example for the v{(j)) dependence, (v — 2 neighborhoods 4, 5 and 8) 
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and for a range like range B (diagonal): 



J 6=6, sin(/) + 



+ COS( 



V^TTtai?^^ + tan 0^+i 

\/l + tan^ (j)i + tan (j)i 



1 i /^^(^ + f) 



V2 tan(^ + f) 



_ 1 / (y^2(l + tan2(/)j+i) - 1 + tan(/>j+i)(l + taiK 
\/2 I (l + tan(/)i+i)(Y/2(l + tan2(/)i) - l + tanc/): 



(9) 



(10) 
(11) 



NB: To get the average speed in that range additionally one would have to normal- 
ize the integrals. 
The average for the whole first octant is the sum 






(12) 



For the deviation integrals 



{v{(t)) - Vavfdct) 



(13) 



the additionally needed integrals are 
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and 
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1.3 Results 
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1 / 1 + tan 0^+1 1 + tan < 



2 \ 1 — tan (j)iJ^i 1 — tan ( 



(15) 
(16) 



However all integrals are simple and analytic, the analytic results do not provide 
too much insight and in the following only the numerical results are given. The 
following average speeds and deviations were calculated: 



Table 1: Average speeds and relative deviations (by angle) for all complete neigh- 
borhoods up to v=10 
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0.79 
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0.105 


4 


1.59 
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...to be continued... 



...continued... 
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However there are ambiguities, these information point to this choice of neighbor- 
hoods (of which only the ones for speeds 1,2,3,5 and 6 can be composed out of the 
corresponding number of subsequent von Neumann or Moore steps): 



speed 


neighborhood (d^„^) 


1 


2 


2 


5 


3 


10 


4 


18 


5 


29 


6 


40 


7 


53 


8 


72 


9 


89 


10 


109 



resulting in this quarter of speed neighborhoods (an agent with maximum speed 
Vjn can reach all cells with a number < Vjn)'- 
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2 A model of pedestrian motion 



As it is not in the main scope of this paper, the model which makes use of the ideas 
above will now be presented only in short: 



Space becomes discretized into quadratic cells with 40 cm as length of an edge. 
Each cell may be occupied by at most one agent. 

Before the beginning of the main part of the simulation, the individual parameters 
- as the maximum speed - are spread over all agents. 

Then the agents are assigned (deterministically or randomly) to their starting posi- 
tion. 

Round by round the agents repeat the following steps until all agents have left the 
scenario via an exit: 

• All agents in parallel choose one of the cells (destination) within the neigh- 
borhood assigned to their maximum speed. 

• The agents sequentially try to reach their destination cell. 

The rules for the selection of a destination cell are quite complex, while the rules 
of movement are rather simple. Former ones are probabilistic. Most important for 
the decision process is the higher probability to select a cell as destination if it lies 
closer to the exit (probability p oc exp(ks{Sjnax — S)), with coupling strength ks 
and distance S to exit). But also herding behavior, inertia, and the distance towards 
other agents as well as walls can play a role. In many aspects this part of the model 
is a higher speeds extension of the model described in Refs. iHEl. 
The rules of actual motion are deterministic, however the sequence in which the 
agents carry out their steps is chosen randomly. Each agent moves within a Moore 
neighborhood to that cell that lies closest to his destination cell. If no cell is avail- 
able that is closer to the destination cell than his current position he remains where 
he is. A once used cell remains blocked for the rest of the round (see Ref. 1 3 1). 

3 Testing the symmetry and the discretization artifacts 

3.1 Walking speeds and travel times of single agents 

To test the benefit of the considerations above for the equality of directions several 
simulations were carried out, where one agent moved a distance of 325 (the number 
< 1000 with the most solutions of Pythagoras: 325^ = A^ + B^) cells into eight 
different directions with two different speeds. Each simulation was carried out 100 
times, ks has been set to 10.0 to make the simulation nearly deterministic. 
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So the overall average evacuation time f or i; = 1 was 302.6 rounds zb 7.00 rounds 
(2.30 %), while f or i; = 5 it was 65.9 rounds zb 0.30 rounds (0.46%). This means 
that standarddeviation of the overall average is roughly by the same factor (5) 
smaller as the speed is larger. If one wants to interpret the agents moving in the two 
examples with 2 m/sec f or i; = 1 one round has to be interpreted as 0.2 seconds 
and f or i; = 5 one round would be one second. Then f or i; = 1 the time to move 
as far as 130 m would vary with the orientation of the discretization axis by more 
than 10 seconds ((326.0-275. 2)*0.2), while for i; = 5 it would be only 2.5 seconds 
(67.4-64.9). 

3.2 A radially moving crowd 

1948 agents were spread over 194812 cells of a circle area (radius 249 cells). With 
four exit cells in the center of the circle the agents started to move at once towards 
the center of the circle. The calculation was done twice: At first all agents had a 
maximum speed 1, during the second run, they had a maximum speed 5. See Fig.s 
[6] for a comparison of how the initially rotationally symmetric spatial distribution 
of agents evolves with time in the two cases. 





Figure 6: Comparison of two simulations with a crowd (black) moving to the center 
of a circle. The left image shows v=l agents after 180, the right one v=5 agents 
after 36 rounds. 



As a last test we compared the simulated walking times of the two alternative routes 

(called A: Start ^2^4^ Exit and B: Start ^ 1 ^ 3 ^ 5 ^ Exit) shown in 

Fig. [71 In reality route B is \/2 times as long as route A and so should the walking 

times be for pedestrians with identical speeds. So this is a specific comparison of 

motion into the two directions and 45 degree. 

These are the average (ten simulations) walking times for agents with a certain 

speed: 
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Figure 7: Two routes: Route A contains horizontal and vertical, route B diagonal 
parts. 
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The deviations from 1 in the last column are due to the integer valuedness of the 
static floor field which leads to probabilistic path-choosing-behavior even in the 
deterministic limt ks ^ oo. For a real valued static floor field the last column 
would contain only 1.00s. 



4 Summary 

In means of minimizing artifacts of discretization, we presented two criteria to 
identify the best neighborhoods for speeds larger one. We presented the results of 
simulations which compared motion in Moore neighborhood steps with motion in 
steps within the best neighborhood f or i; = 5. Depending on the observable the 
results showed the reduction of discretization artifacts by a factor of four or even 
five for the latter neighborhood. This becomes specifically interesting in case of 
finer discretizations, where a subsequent execution of steps within Moore or von 
Neumann neighborhoods would lead to the same dependence of evacuation times 
on the orientation of the axis of discretization, yet on smaller space-scales. 
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